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ABSTRACT 


An  optimization  algorithm  for  a  multi-echelon  model  is  given  which 
does  not  assume  convexity.  The  algorithm  employs  a  bound  technique  to  de¬ 
termine  a  finite  search  area.  Updating  the  bound  is  used  to  increase  the 
efficiency. 
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An  Algorithm  for  Computing  Optimum  Stock  Levels  in  a  Two  Level- 

Maintenance  System 

Introduction 

In  general,  inventory  models  are  optimized  by  some  appropriate  method 
which  depends  upon  convexity  of  the  objective  function.  Thus,  differentia¬ 
tion  and  difference  equation  techniques  are  used  to  find  the  conditions  for 
optimality.  Convexity,  however,  frequently  is  not  seen  where  there  is  more 
than  one  decision  variable.  Nevertheless,  optimization  assuming  convexity 
will  usually  produce  nearly  optimal  values  due  to  the  general  flat  nature 
of  inventory  model  equations  in  the  area  of  the  optimum.  With  this  in 
mind,  an  algorithm  which  assumed  convexity  for  minimizing  the  total  cost  in 
a  two  level  aircraft  maintenance  system  was  used.  During  observation  of  the 
algorithm  it  was  noticed  that  small  changes  in  critical  parameters  some¬ 
times  caused  large  changes  in  the  optimal  stock  distribution.  Even  though 
these  changes  actually  resulted  in  relatively  small  changes  to  total  cost, 
the  inconsistency  of  the  algorithm  was  undesirable  in  view  of  the  time  phased 
nature  of  the  supply  control  study  in  which  it  would  be  employed,  for  over 
the  time  period  some  or  all  of  the  parameters  will  change.  This  paper  pre¬ 
sents  an  algorithm  developed  to  find  the  exact  minimal  cost  distribution 
of  assets  and  eliminate  the  inconsistencies. 

For  a  thorough  discussion  of  the  model  see  reference  1.  An  almost 
identical  model  developed  independently  at  about  the  same  time  as  the  model 
of  reference  1  is  discussed  in  2.  A  sufficient  description  of  the  model 
follows  for  understanding  of  the  optimization  problem. 

Both  models  are  concerned  with  optimizing  stock  levels  of  high  cost  - 
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low  demand  reparable  items  in  a  two  echelon  maintenance  system.  There  is 
a  single  top  echelon  maintenance  facility  which  does  overhaul  and  complex 
repair.  The  lower  echelon  consists  of  several  units  which  are  capable  of 
doing  non  complex  repair.  Demands  for  items  are  placed  upon  the  lower 
echelons.  When  an  item  is  demanded  it;  is  accompanied  by  the  return  of  a 
failed  item.  The  item  can  be  repaired  at  the  lower  echelon  with  a  probability 
f  and  it  can  be  repaired  at  the  top  echelon  with  probability  1-f. 


For  this  class  of  items  the  optimal  policy  at  lower  and  top  echelons 
is  of  the  form  (S,S-1).  Demands  are  assumed  to  be  compound  Poisson  at  each 
lower  echelon  unit.  Under  these  assumptions  an  analytic  formulation  is 
possible. 

The  objective  function  to  be  minimized  is 


where 


N 


N 


TC 


■r 


as.  + 

h  1 


1=0 


I  cbV5i  Wis-V 

i=l 


Sq  «  top  echelon  spares 

Si  =  spares  at  unit  i  of  lower  echelon 

A.f  =  average  failures  per  day  at  unit  i 

^i  =  Percent  total  failures  repaired  at  unit  i 

t'T'i 

Tir=  repair  time  at  i  unit 
Tis=  rePlenishment  time  for  ith  unit 
WQ  =  average  wait  for  stock  at  top  echelon 
=  average  backorders  at  unit  i 
Cjj  =  holding  cost  per  spare  per  year 
Cb  =  backorder  cost  per  spare  per  year 
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TWO  LEVEL  MAINTENANCE  SYSTEM 
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Notice  that 


TC (S  ,S) 
o  “ 


as 

h  o 


N 

i=l 


tchsi 


C,  B.  (S.  \.,f. ,T.  ,W  )] 
b  1  1  i  i  is’  o  J 


or 

TC(So>S)  =  ChSo  +  TCA(So,S) 

where  TCA(S  ,S)  is  total  lower  echelon  cost  given  a  top  echelon  stock 
o  — ' 

of  S  ,  and  lower  level  stock  vector  S  =  (£>.,»  S0,  .  .  .  S  ) 
o  “*  l  z  n 

Any  given  Sq  establishes  the  value  of  the  parameters  of  the  inventory 
level  probability  distribution  which  are  used  to  compute  B^(*)*  With 
parameters  fixed,  B^-)  is  a  convex  function  over  S^.  Thus,  minimization 
of  TCA(So,S_)  can  easily  be  done.  This  is  the  basic  computational  feature 
on  which  the  algorithm  rests. 

Mathematical  Basis  \ 

Before  presenting  the  algorithm  the  following  theorems  need  to  be 
proved . 

Theorem  1.  If  S*  is  the  top  echelon  stock  at  an  optimum,  a  necessary 
condition  for  optimality  is  that  TCA(S*,S_)  be  minimized  with  respect  to 

£  =  (s1,  s2,  .  .  .  sn) 

Proof.  Let  S*  be  the  vector  of  lower  echelon  stock  levels  which 
minimize  TCACS*,^)  and  S/  be  any  other  stock  level  vector.  Now 
TC*,  the  optimum  cost  is 

TC*  =  ChS*  +  TCA(S*,S*)  <  ChS*  +  TCA(S*,S') 

Thus,  S/  cannot  minimize  total  cost. 
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Theorem  2,  Let  S^(j)  denote  the  vector  of  optimum  lower  echelon  stock 

levels  when  S  is  i,  Then 
o  J 

S*(k)  ^  S*(j)  if  j  <  k 

Proof.  Since  each  is  determined  independently  it  is  sufficient 
to  show  that  any  S*(k)  £S*(j),  i-1,2,  •  .  .  N 

When  j  <  k  the  mean  of  the  compound  Poisson  distribution  of  items  in 
the  pipeline  is  greater  for  Sq  =  j  than  for  Sq  =  k.  Let  ji(j)  be  the  mean 
for  Sq  =  j.  Then  jj,  (  j  )  >  jj,(k).  The  condition  for  optimality  at  a  lower 
echelon  unit  is 

G. 

P(SH-1/ jj,)  £  ~  Z  P(S/p,) 

Ub 

where 

CO 

PS/(J.)  =  ^  P(  j  /|J-) 

j=s 

and  pCj/p,)  is  the  compound  Poisson  probability.  Now  P(S/p,(k))  <  P(S/^(j)) 
and  the  optimality  condition  cannot  be  satisfied  with  a  larger  S  if  ^(j)  >  p,(k). 
Theorem  3.  Let  TCA( j ,£*( j ) )  denote  optimum  lower  echelon  cost  when  Sq=  j 

Then  TCA(j,S*(j))  >  TCA(k,S*(k))  if  j  <  k 

Proof.  Let  S*(j)  be  vector  of  optimum  lower  echelon  stock  levels  when 

S  is  j . 
o 

Now  TCA(j,S*(j))  >  TCA(k,S*(j)o  ;>  TCA(k,S*(k)) 

The  first  inequality  follows  from  the  fact  that  backorder  cost  is  reduced 
when  top  echelon  stock  is  increased. 
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Theorem  4.  Let  TC  be  any  total  cost  and  TCA(°o,S*(a>)  be  the  optimum 


lower  echelon  cost  when  top  echelon  stock  is  infinite  (eliminates  W  ) .  Then 
an  upper  bound 

b  =  [  (TC-TCA(ooJS>*(oo)))/c^]  on  top  echelon  stock  can  be 

established  such  that  S*  ^  b. 

o 

Proof,  TC  2s  TC*  -  C^S*  +  TCA(S*,S*(S*) )  >  ChS*  +  TCA(a>,S*(“)) 

Therefore 

[  (TC-TCA(a>,s*(aD)))/Ch]  ^  S* 
where  [x]  denotes  the  largest  integer  ^  x 

A  verbal  proof  of  the  above  is  that  if  you  must  spend  at  least 
TCA(o°,S/fr(<») )  (TCA(«,S*(oo))  is  a  lower  bound  on  lower  echelon  cost)  and  an 
allocation  has  been  found  which  costs  TC,  then  why  spend  any  more  on  upper 
level  inventory  than  TC  -  TCA(oo,S*(oo)  *  There  are  two  features  to  the  com¬ 
putation  of  this  bound  which  enable  recomputation  to  find  a  possible  smaller 
upper  bound.  Suppose  an  initial  upper  bound  b^  has  been  found.  Then  by 
Theorem  2  we  know  the  least  cost  at  the  lower  echelon  will  be  TCA(b^,!S*(b^)) . 
Thus  a  new  bound,  b^  =  [TC  -  TCA(bpj3*(b^) )) /C^]  can  be  computed. 

Likewise,  a  b^  can  be  calculated  as  b^  =  [  (TC-TCAO^  jS^Cbg) ) )  /C^] . 

This  continues  until  there  is  no  change  in  the  bound.  Another  source  of 
improvement  is  from  new  allocations  which  reduce  TC.  Thus  each  time  a  lower 
TC  is  found  the  bound  is  recomputed. 

The  measure  of  efficiency  for  any  optimation  algorithm  is  the  time  it 
takes  to  find  the  optimum.  In  coding  the  algorithm  for  a  computer  we  were 
concerned  therefore  with  time  and  not  necessarily  the  efficiency  of  the 
bound.  The  bound  should  be  reduced  only  if  the  computational  time  in  reducing 
the  bound  results  in  at  least  as  large  a  saving  in  other  computer  time. 
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Quickly,  it  was  learned  that  the  first  means  described  for  reducing  the 
bound  cost  more  than  it  returned.  The  test  algorithm  only  employed  the 
second  method.  Appendix  A  gives  a  listing  of  the  coded  algorithm  which  was 
run  on  a  time  sharing  system  (Corn-Share)  in  the  XTRAN  language.  A  narrative 
description  of  the  algorithm  follows. 


The  Algorithm 


Step  1.  Set  Sq  *  0  and  find  TCA(0,S*(0)  by  incrementing  the  lower 
stocks  independently  until  a  minimum  is  reached.  Set 
TCA(0,S_*(0)  equal  to  current  minimum,  CM,  and  save  stock 


distribution. 

Step  2.  Set  upper  echelon  backorders  to  0  and  find  TCA(ra,S_*(co)) 

Step  3.  Set  bound,  b  =  [ (TCA(0,S*(0) )  -  TCA(^,S*) ) ) 

Step  4.  If  S  s  b  go  to  Step  7 
o 

Step  5.  Increment  SQby  1  (Say  Sq=  j)  and  find  TCA(j ,S*(j)) .  Compute 
TC(j,S*(j))  =  Ch*j  +  TCA(j,S*(j)) 

Step  6.  If  TC(j,S*(j))  <  CM, set  CM  =  TC(j,S*(j)),  save  stock  distribu- 


tion3  recompute 

b  =  [  (TC(  j  ,S*(  j )  )  -  TCA(o°,S*(»)))/Ch],  and  go  to 
Step  4. 

Step  7.  Output  optimal  stock  distribution  and  average  backorders  at 


optimum. 


Step  8.  Stop. 

The  algorithm  must  pass  through  the  minimum  cost  point.  Whereas  if 
the  TC(j,S*(j)>  were  convex  an  optimization  procedure  could  stop  after  ex 
amining  only  one  additional  point  beyond  the  optimum,  this  algorithm  must 
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examine  several  until  the  bound  is  reached.  One  way  of  judging  the  effective- 
ness  of  the  algorithm  is  by  the  number  of  additional  points  it  must  examine 
beyond  the  optimum  before  stopping.  In  testing  of  the  algorithm  over  a  rather 
wide  range  of  conditions  it  was  found  that  the  final  bound  established  was 
always  close  to  10%  greater  than  the  upper  echelon  stock  at  minimum.  This 
was  considered  to  be  good. 

Notice  that  to  insure  reaching  the  exact  optimum,  all  possibilities  of 
upper  echelon  stock  up  to  the  bound  must  be  examined.  Were  the  function 
convex,  more  efficient  methods  of  converging  to  the  optimum  could  be  used. 

A  substantial  amount  of  computing  spped  is  being  lost  to  insure  achieving 
an  exact  optimum,  which  in  turn  insures  a  consistent  solution  pattern  xtfhen 
parameters  are  varied.  Consistency,  however,  was  an  overriding  consideration. 
The  algorithm  is  being  used  in  a  supply  control  study  to  find  the  optimum 
requirements  over  a  5  year  time  span.  Characteristically,  some  parameters 
will  display  slight  changes  over  the  5  years.  The  optimum  requirements 
solution  should  display  a  logically  changing  characteristic  depending  upon 
the  parameter  changes.  If  the  solution  were  to  change  substantially  (this 
can  happen  with  little  change  in  total  cost)  a  good  procurement  schedule 
would  be  impossible.  Moreover,  management  would  have  little  faith  in  the 
solution  if  it  were  to  change  significantly  without  logical  justification. 

This  version  of  the  optimizing  routine  is  not  considered  to  be  final. 
Additional  research  will  be  done  to  make  the  scheme  more  efficient  by  re¬ 
ducing  the  number  of  required  computations. 
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APPENDIX  A 

LISTING  0E  C0UING  E0R  THE  ALG0RITHM 


SUBROUTINE  0PT C UTL*  UL»  SLTL* NB> AMT j  E VAC j AMD AL* C0MUP* ACUP* INV>  W0PT) 
REAL  NSHI  GH..NSL0W 
C0MM0N/ I D/ANAME j DATE 
G0MM0N  /ALL0C/MM,NR»BAGK0 

DIMENSION  AMT  C  N3)  .»  EVAC C NB )  .»  NR  C 2 5 )  .» NR0PT C  25 > 

ANAME=4H0PT 
D  AT  E=  A  HO  9  2  6  6  9 
NSTART=0 
BAVG=UTL 
5  NT0T=O 

BT0T=O.O 
D0  20  1  =  1  >  NB 

AMEAN=AMTC I )+BAVG/UL*EVAC< I )*NSTART 

HACKAR=AMEAN 

c0star=agup*bagkar 

nkai<=o 

10  uagkal=backar 

G0 S  T  AL=  G0  STAR 
NRAR=NRAR+ 1 

OAGKAR=BNAP C  AMEAN i NRAR ) 

I E  C  BACKAR • LE. 0 • ) PAUSE  "ING  AREA  *  BAGKAR=">  #BACKAR 

G0STAR=ACUP*BACKAR+G0MUP*NRAR 

I ECG0STAR.LT. G0STAL) G0  T0  10 

NR  C I ) =NRAR- 1 

BT0T=BT0T+BAGKAL 

20  N TOT  =  NT0T+NR( I ) 

I  E  C NS TART . EQ . 1 ) G0  T0  21 
NSL0W=G0MUP*NT0T+AGUP*BT0T 
NSTART= 1 

G0  T  0  5  , 

21  MM=0 

NSMI  GH=C0MIJP*NT0T  +ACUP*BT0T 
G 0 S T  =  A G UP  * U T 0  T  +  C 0  M  UP  *  N  T0 T 
C0STMN=C0ST 

MMAX= INT  C CNSHI  GH-NSL0  VO/C0MUP ) 

So  MM=MM+1 

BA V  G=  BNAP C  UTL *  MM ) 

NT0T =0 
BT0 T=0 • 0 
U0  50  1  =  1 >  N  B 

AM E AN = AMT  C I )  +  BAVG/UL* EVAG C  I  ) 

NRAR=NKCI) 

BAGKAR  =  BNAP  C AMEAN > NRAR) 

I E  C  BAGKAR . LE • 0 . )P AUSE  SING  NICP  *  BAGKAR= $ > BAGKAR 
CU5TAR=ACUP*BACKAR+C0MIJP*NRAR 
40  NKAK=NRAK-1 

BAGKAL=BAGKAK 
C  0  S  T  A  L = C  0  S  T  A  R 
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LISTING  feJF  C0UING  CONTINUED 


IFCiNJRAR.LT. 0)  G0  T0  45 
BACKAR=  BNAP  C  AMEAN ,  NEAR  ) 
C0STAR=ACUP*8ACKAE+C0MUP*NRAft 
IF CG0STAR.LT. COSTAL) G0  TO  40 
45  NEC  I  >  =NRAR+  1 

NT0T  =  NT0T+lMRC  I  ) 

50  r i  f  0 T  =  BT0  T+  BACKAL 

C03T=ACUP*BT0T+C0MUP*CNT0T+MM) 

I F  C C03T . GE. C03TMN ) G0  T0  60 

C0STMN=C0ST 

NT0PT=NT0T 

M0PT=MM 

B0PT=BT0T 

FOR  1  =  1, NB:  NR0PT C  I  >  =  NRC,I  ) 

MMA X  = I N T  C  C  C0  ST - N SL0 W ) / COM  UP ) 

60  IFCMM.LT.MMAX) G0  T0  30 

FOR  1=1, NB:  NEC  I >=NR0PTC I ) 

<*1M  =  M0PT 

N  i’0  T=NT0PT 

INV=MM+NT0T 

BACK0=B0PT 

W0P  T  =  B  A  C  K  0  /  AM  D  A  L 

RETURN 

END 
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